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Abstract: We present the performance of Multi-gap timing RPCs under irradiation by fully 
stripped relativistic ions (7j8=2.7, Z=l-6). A time resolution of 80 ps at high efficiency has been 
obtained by just using standard 'off the shelf 4-gap timing RPCs from the new HADES ToF wall. 
The resolution worsened to 100 ps for ~ 1 kHz/cm^ proton flux and for ~ 100 Hz/cm^ Carbon 
flux. The chambers were operated at a standard field of £ = 100 kV/cm and showed a high stability 
during the experiment, supporting the fact that RPCs are a convenient choice when accommodating 
a very broad range of ionizing particles is needed. 

The data provides insight in the region of very highly ionizing particles (up to x36 mips) and can be 
used to constrain the existing avalanche and Space-Charge models far from the usual 'mip valley'. 
The implications of these results for the general case of detection based on secondary processes (n, 
7) resulting in highly ionizing particles with characteristic energy distributions will be discussed, 
together with the nature of the time-charge correlation curve. 
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1. Introduction 

The development of large-scale sub- 100 ps resolution RPCs has launched a number of applications 
in Particle and Hadron Physics, notably HARP (PS) and ALICE (LHC), STAR (RHIC), FOPl, 
HADES (SIS 1 8) and CBM (SIS 100) |[T], g g g, |l. The rapid progress of the field ever since its 
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birth [|7|, |8p largely lies on the invention of the multi-gap [|9|], the improvement in the gap definition 
quality that allowed to build an RPC at the 100 ps-a scale for the first time [ |To| ] and the extension 



of the concept to long counters up to 1.6 m []11[]. Detectors of this family are sometimes referred as 
Multi-gap timing RPCs\ hereafter denoted simply by MtRPCs. All the so far existing large walls 
based on MtRPCs emphasize the detector response for mips at moderate particle fluxes (below 1 
kHz/cm^) that are typical environments either in high energy collider experiments or low energy 
fixed-target ones. 

There is, nonetheless, a broad range of newly born applications where MtRPCs need to work 
in highly ionizing environments. Maybe the most remarkable one is the detection of annihilation 
7's for Positron Emission Tomography (PET) [ p^ (based on the detection of the secondary electron 
(by Compton or photo-electric effect)). More recently, in the framework of the R'^B collaboration, 
two main projects started R&D aimed at building large ToF walls for neutron and ion detection, 
respectively [|^, 14 1. Both of them will have to deal with a yet unexplored range of highly ionizing 
particles. 

Remarkable success has been achieved in describing the behavior of these chambers through 
first principle avalanche simulations, revealing the main role of a very strong avalanche Space- 
Charge [15, |T6[ 1^. Despite the presence of this complicated phenomena, a handful of analytical 
expressions can be obtained if standard avalanche evolution is assumed to happen up to the thresh- 
old level [|l^, H^, 19, 1^, 21]. The developed formalisms allow to estimate the influence of the 
ionization loss, hereafter referred in mips units (AE / AEmips) as: 



a,{^E)=K, 
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(1.2) 
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t refers to the average and Oj to the rms of the time response distribution (for the sake of simplicity 
it has been assumed to be Gaussian, so that t can be thought of as the maximum of the time response 
distribution and equals its rms). riti, is the electronics threshold in number of electrons, rig is 
the total number of primary ionizations for mips, g the gap size, a and TJ the multiplication and 
attachment coefficients, the weighting field and Ki is an adimensional constant of order unity 
that contains the effect of avalanche multiplication statistics. The signal rise-time t^ise is related to 
the coefficients of the electron swarm as: 



In 9 



{a-ri)v. 



(1.4) 



in absence of Space-Charge (see [20] for a discussion). The average electron drift velocity is 
denoted by Vg. Typical values of ~ 200-250 ps at threshold level have been measured with 
a careful setup [|^, 23], for typical operating fields of £" = 100 kV/cm, and threshold levels of 



Also simply Multi-gap RPCs or timing RPCs. Nevertheless, it must be mentioned that high time resolution is not a 
property of the Multi-gap technology while high efficiency is also not a property of single gap timing RPCs. 
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<]th = nth<]e — 10 fC. It must be noted that, within the models, eq. LI is generally exact, being eq. 
1.2 full-filed in the limit AE / AE,„ips °° and, contrary, eq. L3 for AE / AE,„jps 1 (otherwise the 
latter must be re-interpreted as a lower limit [ETl]). Despite the development of analytical tools. 



no systematic attempt was done to clarify these dependences with exception of the controverted 
data from [24] and [25], that will be re-visited here, and the much too short survey of [18] (only 2 
points). 

The detection of ions up to A = 200 (Z ~ 100) over large surfaces (~ 5 m^) at relativistic kinetic 
energies Ek = 700 GeV/A has recently brought attention in the R^B experiment, proposed within 
the new Facility for Anti -proton and Ion Research (FAIR) at Darmstadt, Germany [|l^]. In view of 
this potential new application, we conducted systematic measurements at GSI-SIS18 for evaluating 
the detector response under ions up to charge state Z = 6. Complementary, measurements with a 
highly mono-energetic diffuse proton beam (Ek = 1.76 GeV, Oe^/E^ = 4%) and cosmic rays were 
also performed. Spare cells from the HADES system were used [26] and comparisons to previous 
results will be referred when appropriate. 

Importantly, paralell plate geometries (of the kind of PPACs) have been used for many years 
for detecting ions at the 1 00 MeV energy scale or below with great success [ 27 , 28 , 29 ] . Apart from 
extending the time resolution below 100 ps, the introduction of the RPC technology should greatly 
improve the chamber stability and easy its construction, being intrinsically 'spark-protected' at 
ambient pressure. As shown in this work, these features are already provided by typical MtRPC 
designs if ions have enough energy to penetrate in a relatively bulky detector (~ 20%Xo here) and 
ion fluxes are below 1 kHz/cm^. 

The structure of the paper is as follows: the three different experimental setups are explained 
in section ^ the behavior as a function of HV and rate is presented in section ||, while section ^ 
is devoted to the behavior as a function of particle type, a discussion follows in section || about 
the practical use of such a counter in a highly ionizing environment and finally in section ^ we 
summarize our conclusions. 



2. Experimental setup 
2.1 The RPC cells 

Two RPC cells from the new HADES MtRPC wall have been used and allocated inside their cor- 
responding shielding profiles (acting as Faraday cages) in a custom designed aluminum gas box. 
As compared with the HADES system, the gas box was lacking of an optimized PCB for signal 
feed-through of Zg = 50Q. (characteristic impedance) while the routing of the signals inside the box 
was such that they could be read out from the same box side. Electrical isolation of signal cables 
was ensured whenever possible in order to reduce cross-talk. A picture of a typical cell is shown 
in Fig. |l|, consisting of 4 gaps of 0.28 mm thickness, with 2 mm thick aluminum and float glass 
plates. The dimensions of the cells were 22 x 140 mm^ (width x length), corresponding to the low 
polar angle region of the HADES wall. Technical details can be found in [^]. 

The read-out was based on the HADES FEE electronics that provides x40 amplification factor 



in signal amplitude at 2 GHz bandwidth [3C]. The detector was read out in single-ended mode. 



with a signal being taken from the central electrode (anode) after filtering the HV level. No special 
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Figure 1. Up: pictures showing a typical HADES cell, before and after being placed in its corresponding Faraday cage. 
Down: side and front sketch of a HADES cell together with the electrical scheme used for these measurements. 



care was taken for matching the detector impedance (Z^ = 12 it 2^) to the coaxial LEMO cables 
used for signal transmission (Z^ = 50D.), while the FEE built-in dead-time of the order of 50-70 
ns provided stability to the electronic chain by avoiding re-triggers caused by reflections. 2 MBs 



(Mother-Boards) and 2x2 DBs (Daughter-Boards) [31] fed by a customized low-ripple distributed 



LV system [ p2| ] were used for reading-out 8 electronic channels (2x2 RPC cells and 2x2 reference 
scintillators). The RPC -FEE thresholds were set to Vth = 50 mV (equivalently Vth = 50/40 = 1.25 
mV and qfi, ~ 30 fC at the pre-amplifier input) and were not changed during the measurements. 
The gas box was used as a central ground node being the FEE, HV and cathode strips connected to 
it by screws and/or conductive meshes. With this grounding scheme the external noise levels were 
still modest and indeed for Vth < 20 mV some channels were unstable. The digitization and event 



building was done with the Trigger and Readout Board (TRB) developed in [ ]330 and based on the 
HPTDC ASIC chip. 

The oxygen content measured at the output of the gas box was below 200 ppms (measured 
with an Oxygen Transmitter 02X1 of GE Infrastructure Sensing) for a gas flow of 150 cc/min of 
an isobutane-free gas mixture based on C2H2F4/SF6 (90/10). Due to the high electro-negativity 
of the gas mixtures used for timing such ppm levels are not expected to influence the detector 
behavior. 
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2.2 The trigger and reference system 

The arrangement of the RPC cells and reference scintillators is shown in Fig. ^ Two fast Bicron 
scintillators BC-420 and BC-422 read out in both ends by Hamamatsu H6533 photo-multipliers 
(PMs) were used for providing trigger, reference time and track selection for efficiency studies. 
The scintillator-PM assembly was placed inside PVC tubes of 50 mm diameter and optical silicone 
was applied to the contact surface scintillator-PM while the scintillator was wrapped in silver paper 
to improve light transmission. Such a detector choice had previously yielded (Tj, = 35 ps per counter 



for mips [|11[]. 

For the scintillators readout the HADES FEE was also used with the pre-amplifier being by- 
passed in virtue of the PM amplification. The PMs were operated at a nominal voltage V = —2.2 
kV for which the average signal amplitude was just slightly higher than that of RPC signals after 
amplification. Variations around the nominal voltage were introduced to account for slight differ- 
ences in the characteristic gain curves of the different PMs. The same threshold as in the RPC -FEE 
was used. 
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Figure 2. Different views of the arrangement of the reference scintillators and RPC cells used for the measurements. 



We have chosen a setup with the 4 counters (2 RPCs and 2 scintillators) being placed vertically 
along their longest side. The scintillators, being prisms of dimensions 10 x 20 x 100 mm^ (BC- 
422) and 10 X 30 X 120 mm^ (BC-420) were stacked over their narrower faces (10 x 100 mm^ and 
10 X 120 mm^, respectively). They could be aligned with an estimated precision of the order of 
±2 mm. A similar procedure was followed with the RPC cells. Being mechanically identical and 
placed inside aluminum profiles, both RPC cells can be easily stacked with a mutual alignment 
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better than 1 mm. The reference system was carefully centered with respect to the RPCs in order to 
provide a good reference for efficiency estimates. Requesting an additional coincidence with one of 
the ('reference' in such case) RPC is expected to compensate for any error on the aforementioned 



procedure and was done when possible (measurements of section 2.5.2). 

For all the measurements the trigger was provided by a coincidence of the 4 PM output signals. 
In order to do that the PM outputs were split, with one being sent to a Leading Edge Discriminator 
(LED) while the other was fed directly into the FEE after by-passing the amplifier, both discrimi- 
nators having the same threshold. 

Additionally to the scintillator reference system, two fast position sensitive mono-crystalline 
diamonds were placed 14 m upstream, roughly at the focal plane of the last HADES beam-line 
quadrupole. The surface dimensions of the detectors were 4.7 x 4.7 mm^ and 3.5 x 3.5 mm^ 
and their individual time resolutions in the range 100-150 ps [0]. Although not in the trigger 
because of technical reasons, a coincidence probability of 10% was observed off-line, allowing for 
an improved track selection and energy spread determination. Those were used only in a reduced 



set of the measurements (section 2.5.2). 



2.3 Prompt charge determination 

The electron induced (prompt) charge qp was codified in the width of the FEE (LVDS) output 
signal through a 'charge to width' algorithm [pO|], that will be referred as QtoW . The algorithm is 
non-linear as illustrated in Fig. |3| with avalanches lying mainly on a first (steep) linear part of the 
QtoW vs qp correlation curve, while streamers are concentrated in a second (soft) one. Despite the 
non-linearity, avalanches and streamers can be well resolved (Fig. ^-right) and appear separated at 
around <7p ~ 5 pC {QtoW = 200 ns). This 'QtoW method' can indeed accommodate a very large 
dynamic range while keeping the charge resolution below 10%-a for avalanche-like pulses with 
qp > 50 fC. (Fig. |3|-left, dot-dashed line). Points in Fig. ^left were obtained by injecting the 
charge with a fast square pulser (f,-,-^^ = 0.35 ns) after differentiation (C = 1 nF, R = 50Q.) in order 
to emulate the shape of the RPC signals. Data was taken with a Tektronics TDS7104 Oscilloscope 
(BW=lGHz). 

The bi-linear behavior of the average signal width (QtoW) as a function of the average prompt 
charge (qp) can be embodied in a simple parameterization as: 

QtoW = a{l-e-''^")+cqp+ QtoW,„in (2.1) 

after which the resolution of the 'QtoW method' can be obtained: 

(2.2) 



^q,, _ dqp 



qp dQtoW qp 

qp abe^'^^p +c qp 

QtoWmin (53 ns for the channel shown in Fig. ^) is the minimum output signal width, corresponding 
to the minimum time during which the comparator is being self-latched. By numeric integration of 
the signal amplitude, a close agreement for the QtoW vs qp curves was observed when analyzing 
avalanches originated inside the detector as compared with pulser data, while streamers clearly 
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1200 



Figure 3. Left: average behavior of the output signal width (QtoW) as a function of the input pulser charge together 
with a 4-parameter fit. The avalanche and streamer regions are indicated and also the threshold level. The dot-dashed 
line shows the charge resolution resulting from the fit after including the fluctuations in the QtoW value (right axis). 
Right: QtoW distribution measured in the RFC under cosmic rays for an operating voltage V = 6.0 kV. The two separate 
distributions can be attributed to avalanches and streamers. 



deviate from tlie pulser beliavior (up to a factor 1/2 less in the reconstructed qp). The fluctuations 
of the signal width are of the order of Gq,^^^ = 200 ps and no dependence with the input charge was 
observed. A detailed description of the performances of the final HADES -FEE will be the subject 
of a more technical forthcoming publication and will not follow in this paper. 



Eq. ^ provides an accurate phenomenological description and illustrates the different be- 
havior of the algorithm for low and high charges, but cannot be inverted, so in practice a 6''' order 
polynomial was used to obtain the 'calibration curve' qp{QtoW). To avoid errors resulting from 
an incorrect extrapolation, values of the signal width in excess of 290 ns are considered as an 
'overflow' and the maximum pulser charge of qp = 25 pC is assigned to them. Such big charges 
are rarely achieved under ordinary circumstances and they occurred seldom even in the very harsh 
environment studied here. 

2.4 Total charge and rate determination 

The dynamic behavior of an RPC at high rates is bound to the average total avalanche charge q^ 



[ ]35| , |360 . Being readily obtainable by a direct current measurement (after dividing by the avalanche 
rate), its value is directly related to the average gap voltage once the stationary (DC) situation is 
reached, as: 

V = V -IR = V -q^^pd (2.4) 

with <I> the avalanche flux [in Hz/cm^], p the electrode resistivity and d its thickness. V is the 
applied voltage, 7 the average current and R the electrode resistance. As compared with the sta- 



tionary (DC) situation described by eq. 2.4, the interpretation of qj{t) during the stabilization of 
the field in the gap is more intricate [37]. Although it is not the main focus of the present work to 
study the transient behavior of the RPC cells, proper means for determining the average avalanche 
charge as a function of the irradiation time were devised as follows. The particle rate was measured 
from scintillator Si as the coincidence of signals from its two PMs. This provides an unbiased rate 
estimate by suppressing single electron noise. The coincidence signal was sent to the scaler input 
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of a commercial LabJack U3 acquisition board connected to a computer, and its value stored every 
0.2 s. Complementary, the RPC current was measured via the analog output of a 2-channel CAEN 
N471A HV supply, after calibration. Although the HV-display resolution is ±1 nA, by averaging 
the analog output the resolution could be improved down to ~ 0.5 nA. Due to the presence of envi- 
ronmental noise, the measured current had to be averaged over a pretty large time interval of ~ 0.8 
s being stored every 0.2 s, and read-out with the same acquisition board. With the help of a second 
scaler input in the LabJack U3, the cycles of an external 40 kHz clock were counted in order to 
provide a stable time estimate. 

2.5 The physical environment 
2.5.1 C^^beam 

Measurements were performed with three distinct setups. In a first one (Fig. ^a) C'^ was injected 
in the experimental area at a kinetic energy Ek = 1.8 GeV/A directly from SIS 18. The beam was 
focused at the usual HADES target position, 14 m upstream our reference detectors. The latter 
were placed close to the beam dump and optically aligned with respect to the beam-line. A typical 
beam profile at the level of few mm^ was expected in the target transverse (XY) plane, from which 
experience tells that the transverse dimension of the primary beam amounts to some cm^ at the 
beam dump. The accelerator was operated in the so-called slow extraction mode that allows for a 
fairly sustained beam intensity along 8 seconds spill (effectively) with a duty cycle close to 50%. 
Typical spill time-profiles are shown in Fig. |5| (left). 

The non-perfect detector-beam alignment resulted in a position distribution over the reference 
scintillator system characterized by a narrow peak (1 cm— a) close to one of the detectors ends, 
submersed in an uniform background (Fig. |6|-left). After selection based on the signal width in 
the reference scintillators (Fig. the peak could be attributed to the primary C^^ beam while 
the uniform background was largely populated by secondary protons and He together with species 
with Z = 3,4, 5 (either Carbon charge states or Li, Be, B species), the latter at a much lower yield. 
Particles different from Carbon must have been originated along the roughly 14 meters downstream 
from the exit of the vacuum pipe to the detector setup. Due to the different shape of the scintillator 
signals as compared with RPC avalanches, the absence of amplifier and the eventual saturation of 
the PMs at the higher charges, the calibration curve from Fig. |3| could not be used and we present 
in Fig. the raw values of the signal width {QtoW). Based on them we believe, nevertheless, that 
identification can be performed with little ambiguity. 

The position in the scintillators was determined by constraining the width of the time differ- 
ence distribution (left-right) to the known detector dimensions. In such a case the propagation 
velocity can be obtained as: 

Vprop — 2-7- —T (2.5) 
[tL-tR}^Q 

where {ti — t^j^^ refers to the width of the time difference profile at 50% drop from the flat top. The 
effective light propagation velocities obtained were Vprop,Sa = 0.385c and Vprop,si = 0.380c, being 
c the speed of light. Finally, the detectors were mutually aligned by software by imposing that the 
average position difference (y^^ — J.?,) was centered at zero. 
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Figure 4. a) Setup for focused C'^ irradiation (kinetic energy Ek = 1.8 GeV/A, o^^/Ek < 1%). b) Setup for diffuse 
proton irradiation {Efc = 1.76 GeV, o^^/Ek = 4%) after C'^ reactions in a secondary Beryllium target placed 33 m 
upstream the experimental hall. Diamonds have been added and also a second RPC (Ro). c) Setup for cosmic muon 
detection. 
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Figure 5. Illustration of characteristic spill time-profiles. The left column shows the average flux profile (over many 
spills) and the average total charge q^. behavior (right axis) for illumination with C'^ (kinetic energy E/c =1.8 GeV/A) 
at two different beam intensities. Right column shows similar plots for diffuse proton illumination (Ej^ = 1.76 ±4% 
GeV). The extraction time was about a factor 5 smaller in such a case. The transient behavior is little visible since the 
measured current had to be averaged over a pretty large time-interval of 0.8 s due to external noise. 



1 




Y [cm] Y [cm] 

Figure 6. Left: Position distribution along the reference scintillator So under C'^ illumination. The proton, He and 
C'^ curves are shown after charge selection, manifesting the primary C'^ beam (kinetic energy E^ = 1.8 GeV/A) close 
to the border of the reference counters. Right: Position distribution under diffuse proton illumination (E^ = 1.76 ±4% 
GeV) and cosmic ray profile (overlaid). The slightly bigger profile under cosmic rays stems from the larger angular 
spread (see text). 



Only one RPC cell (Ri) was used during these measurements since Rg could not be powered 
up due to technical reasons. A constant voltage V = 5.6 kV was applied to the RPC, while the 
primary flux ranged from 10 to 1000 Hz/cm^. 
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Figure 7. Two-dimensional correlation plot between the raw signal widths obtained in the reference scintillators after 
offset subtraction (color on-line). Six blobs can be identified being the most prominent ones attributable to protons, He 
and C'^. The other 3 blobs with a much lower yield must be Li, Be and B or charge states of the C'^ ion. 



2.5.2 Proton beam 

In a second setup (Fig. C'^ at Ek = 1-8 GeV/A was collided into a Beryllium secondary 
target placed 33 m upstream the experimental area. The magnetic field of the following dipoles 
and quadrupoles in the beam line was adjusted for accepting forward protons from the reaction at 
an average energy Ek = 1.76 GeV (that accounts for energy losses in the Be target itself). As a 
consequence of Fermi motion, protons had still a sizeable angular straggling that rendered a much 
more uniform irradiation than for the C'^ beam (Fig. fright). 

In this setup, auxiliary mono-crystalline diamonds were placed 14 m upstream the scintillator 
reference system. From the time spread between the diamond detector and the first scintillator, 
^T[D-s,) ~ P^' ^^'^ simultaneously measured scintillator and diamond resolutions (a, 
150ps,a,,,^, = a,,,,)=40 
m could be determined as: 



'r(D-5o) ~ f"' ...w . \^T(D) 

150 ps, a^-^^. J = aj-(j, J = 40 ps) the energy spread of the proton beam after propagation over D = 14 



Ek Ek V2 y^l^ tof 
being 



Oek _ mp li^ Otof ^2 6) 



tof=-^ (2.7) 



at„f = Jo^ -a2 -a2 ^ (2.8) 

'"J A/ T{D-Sc,) T(D) T{So) ^ ' 
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resulting in a value for = 4%. The energy losses due to the presence of 14 m of air are much 
below 4% for protons at this energy, supporting that the energy straggling is indeed originated in 
the production at the Be target. Although the time resolution of the diamonds was reduced further 
down to the level of 100 ps after a special analysis [Q], this discussion remains out of the scope of 
the present paper. 

The accelerator was also operated in slow extraction mode but at a reduced extraction time 
(1.5 seconds), resulting in a slightly asymmetric spill time-profile (Fig. |5]-right). Two RPCs (R,,, 
Ri) were powered up and measurements were taken as a function of the particle flux in the range 
10-1000 Hz/cm2 at voltages from 5.0-5.8 kV. 



2.5.3 Cosmic rays 

We took data, after the measurements of the previous two sections were performed, in a 'reference 
setup' with a cosmic stand where the scintillators and RPC cells were assembled according to the 
previous section but horizontal with respect to the ground. From the 2 meters of concrete that cover 
the HADES experimental area it can be estimated that the main source of particles traversing the 
4 detectors are muons at an average energy Ek = 3 GeV (~ IGeV energy loss) [38|. The position 
profile obtained from scintillator So as shown in Fig. fright is fairly uniform but a bit broader than 
for proton irradiation, due to the larger angular straggling (since Si is a bit smaller than So and is 
also in trigger, for perpendicular incidence the measured So length is very close to the one and 
therefore ~ 2 cm smaller than its true length). The HV was varied in the range 5.6-6.0 kV. 



3. Results 



3.1 Performance of the reference system 

The use of LED discriminators for the reference scintillators requires of time-charge slewing cor- 
rections in order to obtain the best timing. The correlation of the time-of-flight between Sg and Si 
(Ats^-Si) "^he average signal width is shown in Fig. ||, being: 



At 



SoSi 



tL + tR 



tL + tR 



So 



(3.1) 



Si 



The different species from Z=l-6 could be cleanly separated as shown in Fig. ^ and a time- width 
correction was performed individually for each specie on the basis of multi-linear segments. After 
this procedure, a combined time resolution of ~ 55 ps was obtained for Ats^-Si that, assuming both 
detectors to perform equally, yielded per counter = 40 ps for C^^, = 39 ps and = 42 ps for 
secondary He and p. The final time-of-flight distributions after corrections are shown in Fig. ^ for 
different particle species, showing a strong Gaussian behavior with ~ 40 ps per counter, while 
the position resolution obtained from the position difference between the counters was ay = 6 mm. 
No dependence with the particle specie was observed. Being the values for Oj consistent with the 



electronic resolution of the mean-time obtained with pulsed signals [ |32| ] and not depending on the 
primary ionization, provides a strong evidence that the performance of the reference system was 
indeed limited by the TDC resolution. 
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Figure 8. Up: multi-linear fit used for correcting the walk of the time of flight between the reference scintillators 
(Ats^-Si) ^ function of the signal width (QtoW) for p. He and C'^. Down: distribution of times of flight between 
scintillators after the correction procedure for the three species. 



3.2 RPC performance 

The track selection was based on graphical cuts on the 2-dimensional signal width distribution 
in the scintillators (Fig. ^ for selecting the particle species corresponding to Z=l-6. Due to the 
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energy regime, either secondary forward protons from C up-stream reactions or primary ones 
with Ek — 1.76 GeV can be highly regarded as mips (yjS ~ 2.7), the latter having additionally 
very little energy spread. On the other hand, an estimate of the energy loss from the Bethe-Bloch 



formula (eq. |4TT1 ) indicates that cosmic muons with Ek — 3 GeV (/jS ~ 30) can ionize 30% more in 
the RPC gas than mips. A normalized charge distribution after cuts is shown in Fig. ^ where indeed 
a very small difference is seen between the 3 different cases, being the average charge released by 
cosmic muons a bit higher, indeed, than that for mips. 
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Figure 9. Comparison of the cliarge distribution measured in the RPC under a diffuse proton beam at a kinetic energy 
Ex = 1 .76 GeV ±4% (continuous line), secondary forward protons from C'^ interactions at Efc ~ 1 .8 GeV (dotted) and 
cosmic muons at Ex 3 GeV (dashed). 



After particle identification, the last 1 cm from the border of the scintillators was disregarded. 
Furthermore, a 'co-linearity' cut of ±2-0 in the distribution of the position differences between 
scintillators (j.,^ —ysi) was applied for improving the track quality selection. In order to reduce the 
effect of the non-uniform irradiation, C^^ ions were additionally selected by a cut at j > 2 cm and 
other species by j < cm. After these cuts for enhancing the purity and quality of the reference 
tracks, the time of flight between the two overlapping RPCs {AtR^^R^) was defined: 



tL + tR 



tL + tR 



Ro 



(3.2) 



and the time resolution per counter was obtained as = cTa/^_^_^^ /V2. Correction curves with 
respect to the charge of both RPCs (slewing correction) and the position as given by the scintillators 
were performed. This is discussed in detail in section A typical time-of-flight distribution for 
protons at an uniform flux of 150 Hz/cm^ is shown in Fig. IT^ yielding = 87 ps per RPC. We 
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applied the usual recursive fit around itl.S-a to reduce the influence of tails in the fit [II], but for 
this set of measurements the value of changed little indeed if a direct fit was performed. 
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Figure 10. Time of flight distribution between RPCs for a diffuse proton beam (kinetic energy Ek = 1.76 GeV ±4%) 
at 150 Hz/cm^ and V=5.6 kV. The inset shows the same distribution in logarithmic scale, indicating the presence of tails 
at 1% level. The counter resolution from the fit is Gj = 123/v^ — 87ps. 



In the C^^ experiment, due to the absence of the RPC cell Rg, the time difference with respect 
to scintillator So was calculated in a similar manner: 



tL + tR 



tL + tR 



(3.3) 



So 



and the RPC resolution obtained after subtracting the contribution of the reference system: 
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3.2.1 Charge distribution 

Normalized charge distributions are shown in Fig. |n| for Z = 1-6 at V = 5.6 kV at low rates (0 < 10 
Hz/cm^). A phenomenological fit to 4-parameter Landau-like distributions 



: aexp{—bqp — ce 



-dq,, 



(3.4) 



provides a good description and partially captures the high charge tails. Fig. [11| shows that there 
is indeed a strong correlation between the charge initially released and the final avalanche charge. 
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Figure 11. Up-left: RPC charge distributions for fully stripped ions with atomic charges Z = 1,2 together with a fit to 
a Landau-like function. Up-right: Z = 3,4. Down-left: Z = 5,6. Down-right: Fitted curves for Z = 1-6. The normalized 
momenta is yji ~ 2.7 in all cases, V=5.6 kV and <p < 10 Hz/cm^. 



but clearly not proportional, as one would expect if Space-Charge would be strong. Nevertheless, 
specific experimental situations where protons must be separated from Carbon, for instance, can 
be certainly addressed with these counters. These issues are further discussed in section ^. Re- 
markably, it seems not to exist a clear separation between avalanches and streamers at high initial 
ionizations (contrary to the case when increasing HV, see Fig. ||). This fact cannot be firmly at- 
tributed to the RPC dynamics since our charge calibration is sensitive to the signal shape and is 



therefore different for streamers, as pointed out in section Due to that, the region > 5 pC 
can be probably not assessed precisely in this work. 

3.2.2 Dependence with HV 

A scan in HV was performed for Ek = 1.76 GeV diffuse protons and cosmic muon irradiation. A 
particle flux of 150 Hz/cm^ was chosen for the former while a much lower one was observed for 



the latter. The behavior of the time resolution and efficiency is shown in Fig. 12. Closed an open 
circles represent the efficiency when at least one or when both detector ends collect a valid signal, 
respectively, showing little difference. As compared with [ |2^ ] the behavior is slightly (but consis- 
tently) worst. While the resolution at the plateau is a^. ~ 90 ps both for protons (diamonds) and 
cosmic muons (triangles), and the efficiency e = 92% (squares) and e = 96% (circles) correspond- 
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Figure 12. Efficiency for cosmic muons as a function of High Voltage when requiring a valid hit in least one detector 
end (full circles) or in both (open circles), together with the efficiency for = 1 .76 ± 4% GeV protons (squares). Time 
resolution for cosmic muons (triangles) and Ex = 1.76 ±4% GeV protons (diamonds) is also shown with scale on the 
right axis. A comparison with the system values from ^ Q (efficiency, resolution and working point) is indicated by 
dashed lines. 



ingly, the values reported for the HADES in-beam test [ Pq ] were a^. = 73 ps, e = 99% (intrinsic) 
and £ = 98% (system). Nevertheless, in such a case, no particle identification was available and ac- 
tually a fraction of the impinging particles is expected to be more ionizing than mips. This, together 
with the slight redundancy of the HADES system (~ 30% cell overlap) may explain the slightly 
worse efficiencies observed here for mips. The worse resolution may be attributed to the grounding 
scheme and overall noise, routing and signal feed-through (unproperly matched) as compared with 



[26]. The multi-peak structure observed in the time response for the lowest voltage run {V = 5.0 
kV, £=24%) suggests that the noise level was abnormally high, eventually influencing the response 
when the signals were very close to the threshold. This effect disappears at the nominal voltages 
V = 5.6 — 6.0 kV but a residual contribution to the overall resolution can certainly not be excluded. 

The chambers showed dark rates of 0.05 Hz/cm^ {V = 5.6 kV) and 0. 15 Hz/cm^ {V = 6.0 kV). 
Few minutes were needed after applying the HV before such values were reached, exceeding those 
by a factor 10-20 otherwise. Dark current was generally at the level of 1 nA or below. 

3.2.3 Dependence with particle rate 

The rate capability of 4-gap MtRPCs with metallic electrodes and 2 mm glass plates has been 
studied under mip irradiation before ||40|], yielding a 5% efficiency drop (— A£) at 350 Hz/cm^ and 
slight deterioration of timing performances. Another study in [pT|],[|42|] consistently showed that a 
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moderate 10 °C temperature increase allows to keep the time resolution up to 1 kHz/cm^ at least. 
We present once again the behavior of the time resolution and efficiency as a function of the particle 
flux in Fig. |l^ for a typical field E = 100 kV/cm, together with the fit from previous work [ |40| ] 
(dashed lines) and added the behavior under C'^ ions. In this last case the non-uniform irradiation 
is taken into account from the measured position profile (Fig. ^ through an average over it 1 cm 
around the peak value. For the proton case the study was performed with a diffuse beam (section 



2.5.2) 



Indeed the behavior for protons as seen in Fig. 13 is similar to the one reported in [40] but 
with a slightly worst resolution at low rates. Clearly, the efficiency for C^^ is unaffected even 
at 600 Hz/cm^ since the effective field drop is overcompensated by the largest initial charge (see 
Fig. ^). In this particular situation the resolution is expected to be more sensitive to the effective 
field (eq. |1.2| ) through trise{E) and so a deterioration of the rate capability for C^^ by a factor x5 
with respect to protons can be inferred from the figure (if the flux at a,- = 100 ps is compared, for 



instance). Note that from eq. one expects (at moderate rates, when q.^ changes little) this factor 
to be directly the ratio of the total avalanche charges q^^^il^Tp- If '^^e MtRPCs would work in the 
proportional regime this ratio would be 7} = 36. 
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Figure 13. Time resolution Oj and decrease in efficiency — Ae in tiie RPCs for protons at = 1 J6±4% GeV (circles) 
and fully stripped Carbon ions at = 1.8 GeV/A (squares) as a function of their flux. The operating voltage was V = 
5.6 kV. 



A remarkable fact is that for MtRPCs under mip irradiation both e and begin to deteriorate 
at a similar value of the flux while for higher initial ionizations the former can be kept at reasonable 
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values even at high particle fluxes and mainly a^. is affected. 

Fig. [TJ shows that the detection of relativistic ions up to Z = 6 is feasible with 4-gap MtRPCs 
operated at typical fields {E = 100 kV/cm), with a resolution ~ 80ps, and e ~ 100%. A more 
detailed differential study as a function of the energy loss is presented in the next section. 



4. Dependence with particle type 



4.1 Introduction 



At fixed momentum per nucleon (yjS constant), heavy charged particles lose energy due to elec- 
tromagnetic interactions in square proportion to its electric charge according to the Bethe-Bloch 
formula: 



dE / Z„, \ 7 1 



1 2m^£2^Y7^_n2 

2 /2 



(4.1) 



where Z,„, A„, are the atomic and mass number of the medium, nie is the electron mass, / the mean 
excitation energy, and K = 0.307075 MeVg^'cm^. A prescription for calculating Z,j,/Am and / in 
mixtures according to [38, ^] will be used in the following. 



Using fully stripped ions of charge Z is a practical alternative to study the energy loss de- 
pendence of a counter over a very broad dynamic range. Moreover, as long as measurements are 
performed in the forward beam direction any low-Z fragment is likely to proceed from spallation 
reactions, roughly keeping its energy per nucleon (Ek/A) and so j8. Therefore in the present exper- 
imental situation, having a primary C^^ beam with Z = 6, it is reasonable to assume that ions with 
Z=l-6 will have a primary ionization in a relative proportion close to 1, 4, 9, 16, 25, 36. Addition- 
ally, for the energies Ek/A = 1.8 GeV used here, even in the very unlikely case that a secondary 
proton (for instance) would travel forward with a kinetic energy as low as Ek,p = 1/4{Ek/A), it 



would release only 20% more energy than mips according to eq. 4.1. The above considerations 
are true only for electromagnetic interactions with hadronic interactions expected to happen only 
residually. In our setup we have no way to separate contributions from hadronic interactions but 
with another reference detector placed RPC downstream, one could separate those by vetoing on it. 
We will disregard the effect of hadronic interactions for the latter discussion and so, being yjS ~ 2.7 
it will be useful to re-express the Z-dependence as AE/AE,„ips = Z^. 

4.2 Charge distribution 

To the authors knowledge, there is only one published work on the energy loss dependence of 



MtRPCs [24], whose conclusions yet remain disputed, specially regarding the very high shifts 



observed for the average time of flight t as compared with the expected one [)25|]. We revisited 
the situation starting from the dependence of the average prompt charge cjp with the energy loss. 
For that we took the average charge for protons as a function of the momentum p from Fig. 8 



in ref . [ 24 ] and pad-ring 3 that corresponds to perpendicular incidence, according to the authors. 



The Bethe-Bloch formula for the HARP gas mixture was evaluated from eq. ^^Jj and used for re- 
calculating AE /Ax, re-normalizing its value to the value for mips. Since the average prompt charge 
is referred in arbitrary units in [^] we arbitrarily re-scaled it to provide the best agreement with 
our data. Note that the geometry of the HARP cells is very close to the HADES one, having also 4 
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gaps and a very similar gap spacing (g = 0.3mm and g = 0.28mm, respectively). Values for HARP 
are shown in Fig. 14 -left (squares). 




Figure 14. Left: average prompt charge qp (circles) and charge at maximum ^p(max) (triangles) from the RPC as a 
function of the energy loss expressed in mips units. Data from HARP is shown by squares for comparison. Different 
trend curves are indicated. Right: average prompt qp and total qj. charge (circles, squares) as a function of the energy 
loss expressed in mip units. A linear fit is shown (dashed line) together with the same curve re-scaled by a factor 6.5 to 
meet the qj. data. 



When increasing the energy loss the probability of streamers also increases [^]. The above 
fact makes the interpretation of the charge spectra difficult since a theoretical description of stream- 
ers in RPCs is not yet available. In order to simplify this task we used the fits to Landau-like 



functions of Fig. 11 to extract the charge at maximum ^^(max) whose energy loss dependence is 



expectedly dominated by the avalanche dynamics and little dependent on the streamer mechanism. 



The values for ^p(max) together with the direct charge average qp are shown in Fig. 14-left. The 
trends of ^^(max) and qp are indeed similar up to AE/AE,„ips = 4 but diverge afterwards, pointing 
to the influence of streamers for AE/AE,„ips > 4. Note the reasonable good agreement in the trend 
of (jp when comparing with the HARP data. The following phenomenological curves could be 
found: 



/ AE 

^p(max)oc J— (4.2) 

qp=qp{\)+0.2{-^) (4.3) 

As a result of Space-Charge, neither qp (max) nor qp show a proportional behavior with the primary 
ionization. 

Since the total current and particle rates were continuously monitored (Fig. ^ the total average 
charge q^ could be determined to be (7j, = 5 ± 1 pC under proton irradiation and q^ = 53 ± 10 pC 
for C^^ under C'^ irradiation. The latter calculation requires to re-estimate the proportion of C'^ 
ions in the measured rate, and is done from Fig. H, yielding 3/5 of the total. The contribution of He 
and proton ions to the measured current is neglected. Note that if we assume the same scaling with 
the energy loss for q^ than for qp the He and p rates (2/3 of the total) contribute to the measured 
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current in less than 10%. From run to run variations we estimated in 20% the error in for protons 
and the same uncertainty was taken for the C^^ value. Data for qp and q^ is shown as a function 
of t}£ I t}£mips in Fig. |l^-right together with a linear fit to the qp data. The same function can 



describe the trend of the q^ points after multiplying by a factor x6.5. If a model of the kind of [ 15] 
would be used for describing avalanche multiplication under strong Space-Charge, the ratio qjlq^, 
is indeed expected to depend on the initial energy loss, because the higher the initial charge the 
earlier the avalanche enters in the Space-Charge regime and the more the released electrons drift 
before arriving to the anode, resulting in a higher electron-induced charge for the same avalanche 
charge. Because of that, and because of the different behavior of qp and ^p(max) with AE/A^m/ps, 
we believe that the scaling observed in Fig. |I^-right is strongly affected by the presence of the 
streamers. 

4.3 Time distribution 

A striking observation in [|l|] was the existence of very strong drifts of the average time-of-flight as 
a function of the energy loss. The much larger range of ionizations present here allows us to discuss 



this point in high detail. First, in Fig. 15 -left the time resolution after slewing correction is shown 



for the different species at V = 5.6 kV. The analytical scaling from eq. 1.2 is illustrated by fixing 
OjiX) to be 80 ps. Unfortunately, the scaling of is not visible in data, pointing to the fact that a 
Umitation to the resolution different from avalanche dynamics is present in the measurements. In 



connection with this observation it must be recalled that despite the solid derivation of eq. 1 .2 there 



is no systematic experimental study so far on the behavior of with the energy loss. Proving (or 



disproving) the scaling of eq. [L2| by systematic measurements would be extremely important in 
order to better understand the practical limitations of these counters for timing. An experimental 
confirmation of this scaling would allow to understand why, in practice, the difference between 
1-gap counters and multi-gap can be as little as Oj = 55 ps (1 gap, [p3|]) and = 40 ps (10 gaps, 
[0]), when operated in very similar conditions of gas mixture and field. No data from the HARP 
collaboration is available on this important aspect. Since our system is limited to a resolution of 
= 40 ps for the mean-time (Fig. ^ a much more accurate FEEh-TDC chain would be needed, 
nevertheless, in order to evaluate the theoretical prediction for C'^ ions. 



Fig. |15|-right shows the average behavior of <AtRPc-s„> (e<l- p.3| ) before corrections and 
referred to the case of mips (circles), together with the analytical scalings from eq. LI assum- 
ing t,-ise = 300 ps (dashed line). Despite the reasonable value obtained for t,-ise this is indeed the 
minimum value consistent with the observations: if a sizeable drift is present in the reference 
scintillators part of the drift originated in the RPC would be effectively shadowed after taking the 
difference. In the absence of constant fraction discriminators the best practical approach to es- 
timate this effect in the scintillators is by looking at the Ats„-s, distributions when the particle 
releases charges corresponding to Z=l-6 in So and keeping a narrow cut around Z=l in Si. The 
above value can be used to correct <Atjipc-So > froni the scintillator drifts and provide an estimate 
of the time drift caused by the RPC alone (triangles). The latest HARP data from [ p5[ ] are also 
overlaid (squares) after an arbitrary shift. According to the authors the values have been obtained 
after a re-analysis based on a direct determination of the momentum p by physical constraints 
(elastic scattering) and are expected to be free from artificial drifts that previously arose from an 
incorrect momentum reconstruction. In order to compare with present data the energy loss was 
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Figure 15. Left: Time resolution CT^ as a function of the energy loss expressed in mips units together with the approx- 
imate 1 / \/AE scaling. Right: average time of flight between RPC and scintillator So as a function of the energy loss 
expressed in mips units together with the approximate ln(l/Aii) scaling before (circles) and after (triangles) correcting 
for the scintillator drift. HARP data is shown as full squares. Two theoretical curves are fitted in each case with a single 
parameter frise = 300 ps (dashed line) and trise = 700 ps (continuous line). 



re-calculated from eq. 4.1, as in previous sub-section. Several things can be noted: i) despite the 
different FEE band-widths the corrected values for AtupcSo agree with latest HARP data and can 
be described with a logarithmic dependence if t^ise = 700 ps; ii) surprisingly, the authors in [25] 
implicitly assumed a value for trise = 100 ps for the discussion in the paper, presumably inspired 
by theoretical values for (a — T])vg, despite this seems to be largely inconsistent with their own 
measurements; iii) the previously reported '500 ps' effect from [24] can also not be confirmed by 
the present measurements; iv) a lower limit for trine > 300 ps readily arises from the present data. 

For the sake of completeness, fig. |l^-right shows the detector efficiency as a function of 
the energy loss. The fact that for Z > 1 the efficiency is almost 100% indicates that the trigger 
geometry is reasonable, being the lower values for Z = 1 partly attributable to the slightly low 
operating voltage (Fig. [T^. Fig. |^left presents the specie population after cuts. 



5. Discussion 

5.1 Element identification 

When aiming at ion identification a time of flight measurement provides the ratio = p/A while 
momentum (p) and element (Z) identification must be done by independent means [|T^]. So it is 
interesting to see at which extent a MtRPC has resolving power in Z by using the qp information. 
In order to evaluate this we define the purity P against ion Z2 for a 90% identification efficiency of 
ion Zi by: 

P^^(^] (5.) 



where the condition £zi = 90% imposes a fixed cut in the prompt charge. Taking all the populations 



from the normalized distributions of Fig. 1 1 the suppression factor ITgo for specie Z2 when aiming 
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Figure 16. Left: population of ions with charge states 1—6 after cuts ('good events'). Right: detection efficiency as a 
function of the energy loss expressed in mips units. The points have being obtained for V = 5.6 kV at low particle rates 
((j) < 10 Hz/cm^). 



at detecting Zi with 90% efficiency can be defined as: 

'Nz^ 1 

rigo is shown in Fig. |l^ for the different ions. 



n 



90 



(5.2) 




Figure 17. Suppression factor of element in charge state Z2 when imposing a cut in the RFC charge for a 90% 
identification efficiency of element in charge state Zi . The suppression factor has been obtained after fitting the charge 
distributions to Landau-like functions. 



The overall suppression is modest, but a strong suppression already exists for Z = 4 ions 
when aiming at proton detection (~ x 65). Contrary, the suppression of alpha particles for proton 
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detection as well as that of Z = 6 ions for Z = 4 ion detection is only slightly higher than 1. 



5.2 Time-Charge correlation 



Slewing correction is the procedure after which the time of flight walk as a function of the signal 
amplitude is corrected for. The nature of this correlation is not yet clear, but it seems to be domi- 



nated by the electronic response at low charges []16[]. Indeed very little correlation, if any, remains 
in the present case for Z = 6, where the signals are largely above threshold (Fig. [l8|-up), unlike 
protons. Moreover, the measured time for high charges is systematically shifted by roughly 400 
ps when comparing protons and Carbon, as already discussed in connection with Fig. |l^. In this 
particular case the external detectors provide a clean separation in Z so that the time distributions 
corresponding to the different ions can be identified and the drift of the average time becomes a 
trivial systematic effect. Now, if Z is not known (or, generally, if t}£ j ^mips is not known), this 
systematic dependence of the average time can not be fully corrected for and these variations enter 
effectively in the measured response function. 

The worsening of the time resolution due to variations in the average time due to different 



primary ionizations is known since time [ |19| ] but has only been applied to account for the Poisson 
fluctuations in the initial number of clusters. Following the spirit of this work, a simple although 
more experimentally-driven approach is devised here to estimate this effect in practice. Let's as- 
sume first that the time resolution does not depend or depends little on the average ionization. This 
assumption is not supported by avalanche models, but would be the case if the measured reso- 
lution is not dominated by the detector response function but by an external source (FEE jitter, 
noise and/or TDC resolution). We think that this case is practically relevant. Further, we assume 



a flat distribution of possible energy losses over a range AE,, 
distribution can be estimated from the relations: 
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(5.3) 

(5.4) 
(5.5) 

(5.6) 



mips 



Eq. 5.6 has the curious property of tending asymptotically to f,-,-ig/ln9 for high energy deposition. 



Since the intrinsic detector resolution is expected to be of this order (eq. 1.2) in the limit of low 
ionization, one may conclude that energy spread cannot modify the detector resolution by a large 



factor. In reality two things can happen: i) that the intrinsic detector resolution from eq. 1.2 is 
reached, and then from eq. ^j6| the worsening due to the energy loss will be of the same order 
than the resolution itself, being the latter partially compensated by the improvement in the intrinsic 
resolution for every value of AE (eq. |1.2|) or ii) that it is not reached, and then the deterioration 



from eq. |5^ will be smaller than the resolution itself, anyway. So it seems that the variations on the 
initial ioniztion are called to be 'second order' for the counter resolution. The above considerations 
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Figure 18. Up: Scatter plot (color on-line) showing the walk of the time-of-flight between RPC and scintillator So as 
a function of the signal width (QtoW) for protons (o), Helium (x) and Carbon (•). Tendency curves are shown in each 
case. Down: time-of-flight distribution after applying an unique slewing correction for all species (ignoring external 
identification in Z). Identification is applied at the end of the procedure to evaluate the result in such conditions. A 
resolution CT^ = 1 15 ps is obtained for the combined fit, with increased tails towards delayed times as compared with the 
case where identification in Z is provided. 



cannot account for all the possible AE" distributions and must be taken with care when a highly 
Gaussian response is needed. Relevant cases where this effect should be taken into account are 
neutron or 7 detection due to interaction in the electrodes, where a secondary particle with its 
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corresponding energy distribution is released, and also ion detection. So, the usual assumption that 
a higher ionization renders a better time resolution can not be considered universally true for this 
kind of detectors, on the basis of the existing knowledge. A detailed discussion on this effect can 



be found in [|44|]. 

Another relevant situation would be the dependence of the average ionization with Z, assum- 
ing the energy per nucleon to be the same for all the particle species. Let's further assume that 
2N + 1 species are present in the same amount, centered around specie Z. If a detector providing 
identification in Z would not exist an extra jitter will appear in the form: 

W l^\n{Z + if + In (Z - if ^ InZ^ \ 




lN+\ lN+\ 
^2N (ln(Z + /)2-Mf )'+('in(Z-/)2-M^ 

- []^) L ^ (5-8) 



in the hmit where Z and A'^ are big, but still N,eq. 5.8 yields 



so the practical influence of the uncertainty in Z in the measured resolution is also little for high 
Z values. However, for cases where Z ~ A/^ ~ 1 (the present case) the jitter arising from a bad 



identification in Z will be as high as 500 ps, as illustrated in Fig. 18. As a benefit, the Z-resolving 



power of the RPC itself will be also much higher, resulting in a time-charge correlation curve 



different from the usual one for mips (compare circles in Fig. |18| with all the three species together). 
In particular, if only Carbon and protons would be present in our sample they can very easily 
identified and corrected for, and only few events (Carbon with QtoW < 100 ns and protons with 
QtoW > 100 ns) would be wrongly identified. If we would have ignored for a while, in this 
particular experimental conditions, the existence of secondary particles arising from the primary C 
beam, and would have used a single correlation curve the net result would have been a deterioration 
of the estimated time resolution for C ions from a^. = 73 ps to a^. = 1 15 ps, with large tails towards 
delayed times. 

5.3 Space-Charge 

There is a simple analytical model that accounts for Space-Charge effects, and that is based on the 
avalanche equations in the presence of a charge-dependent Townsend coefficient |[l^]. This model 
is analytical for each avalanche and depends on its position and initial ionization, but the solution 
to the general problem can only be treated numerically. It is not the purpose of the present work to 
quantitatively describe Space-Charge, but we did some clarifying attempt by evaluating formulas 



in [15] for three typical avalanche positions in the gap. Further, the simplifying assumption that 
all the charge is released at the position of the avalanche is done. We took values a* = 73 mm^^ 
?io/gap= 3AE/AEmips, qsat = 10^ (number of electrons for a 50% drop of a*) and leave the overall 
normalization as a free parameter to match the data. Finally, we looked at three situations (Fig. 



19): i) the avalanche is produced in the cathode {x„ = 0, dotted line), ii) in the center of the gap 



{Xo = g/2, dashed) and iii) 20/xm separated from the anode (dot-dashed). 
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Figure 19. Measured charge at maximum ^^(max) as a function of the energy lost in mips units (circles). The \/AE 
experimental scaling is shown together with the evaluation of the model of Fonte [ [1^ ] under the assumptions of the text, 
for three different positions Xo of the avalanche inside the gap (center, close to the cathode and close to the anode). 



Clearly, we cannot conclude from the comparison shown in Fig. |T9|by such a simplified model, 
but we can illustrate how Space-Charge works. For avalanches produced close to the cathode there 
is almost no sensitivity to the initial ionization, since its evolution is anyway doomed to a similar 
end once the critical value of qsat is reached. For avalanches produced close to the anode, the 
resulting charge will be obviously rather small (an effect absorbed here in the overall normalization) 
but essentially proportional to the initial one, since the avalanche charge is then much smaller than 
qsat and Space-Charge plays no role. For the center of the gap, an intermediate situation is reached, 
where the avalanche growth is modified but some sizeable correlation with the initial charge is still 
present. Remarkably, this latter situation matches well the trend of the experimental points. 

5.4 Streamers 

From the experimental point of view when aiming at ion detection it would be probably convenient 
to operate the detector in absence of streamers. It is possible to perform a semi-quantitative evalua- 
tion of what is the practical implication of this fact by using some classical arguments: if we denote 
by a* the effective gain threshold at which streamers would arise, then following the approximate 



derivation of [ |450 it can be seen that this value will depend on the initial ionization as: 

InZ^ 

<^c\z = (^c\ mips — (5.10) 

SO, the critical value of the Townsend coefficient for an ion with charge Z is naturally lower than 
that of mips, and the same will happen for the field at which streamers will appear in one case and 
the other. Now we can further assume a dependence of a*{E) ~ — b (inspired by transport 
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codes [|lq]) and obtain the necessary decrease in the field to operate the chamber in the absence of 
streamers that would simply be, for a given ion Z: 



£[kV/cm] = ^ [a:Up, - l/glnZ' + b] 
if we make the natural assumption a*\mips ^ cc*{Eo = 100 kV/cm) then: 



(5.11) 



dE 

^^[kV/cm] = £o[kV/cm] - — 



lnZ2 



(5.12) 



as a function of Z. By directly substituting the simulated value for [16] it can be expected that 
even for the most extreme case of Au detection (Z = 79) a mere decrease of a 10% in the field 
would be sufficient to operate the chamber in pure avalanche mode, and the anticipated decrease 
of performances resulting from the lower field would be over-compensated if the scalings with Z 
illustrated through eqs. L2 and 13 hold: 

' Z ln9 



£(Z) > 1-exp 



1--!- 
a 



ag 



(5.13) 
(5.14) 



since t^ise and a have approximate linear dependences with the field E. 



6. Conclusions 

A time resolution Oj ~ 80 ps at e ~100% has been consistently measured for relativistic ions with 
charge states Z=2-6 for the first time by just using standard 'off the shell' Multigap timing RPCs 
from the HADES wall. 

The energy loss dependence of the avalanche charge and detector-related time-of-flight sys- 
tematic shifts have been compared with previous data and extended over a much larger range of 
primary ionizations, showing a reasonable agreement. The measured time drifts cannot be accomo- 
dated in the existing theoretical framework unless a value for the signal rise-time 2.5 times bigger 
than current experimental estimates is assumed. The observed behavior of the prompt charge with 
the initial ionization seems to provide a stringent benchmark for the parameters of Space-Charge 
models and a simple comparison was attempted, showing a reasonable agreement. 

Operation up to Carbon fluxes of 100 Hz/cm^ was demonstrated with < 100 ps and e ~100% 
under an 8-second spill (50% duty cycle) and ~ 2 x 2 cm^ irradiated area. Above 100 Hz/cm^ the 
time resolution deteriorates rapidly but the efficiency was kept up to 600 Hz/cm^ (at least) due to 
the much higher initial ionization as compared with mips. The behavior of the resolution as a func- 
tion of the flux is similar to that of mips, when re-scaling the rate by a factor x5. This value is very 
similar to the measured ratio of the total charges q^^^^J = 6.5, as expected from a simple DC 
modelling of the detector. Based on this observation, the measured trend of q^ (Z) can be extrap- 
olated to high Z, yielding an approximate dependence for the rate capability as ^,„ax{Z) = (|^? , 
that would severely limit the operating rate to barely 1 Hz/cm^ for Au ions. Based on the present 
measurements, it is likely that the very high initial charge will not affect the chamber stability in 
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such situation but it will largely reduce its rate capability. In a realistic application, the working 
voltage should be chosen taking this fact into account. 

The nature of the time-charge correlation and the practical timing limitations when the energy 
loss of the ionizing particle can not be addressed by external means have been discussed. It has been 
shown that under reasonable assumptions the detection of particles through secondary processes 
(in case of neutron or y-photons, for instance) would yield an extra time jitter of the order of the 
detector intrinsic resolution for mips, therefore the good timing characteristics of these devices 
will be preserved even in such a situation. A practical example has being given based on the 
present data. Nevertheless, the intuition that higher ionization yields better results does not seem 
to be a trivial statement for these counters and every physical case should be probably addressed 
experimentally. 

Despite the large dynamic range explored here the chamber+electronics performed stably dur- 
ing the whole experiment, underlining the superior performance of float glass Multi-gap timing 
RPCs in highly ionizing environments when high rates are not required. Always depending on the 
Z of the species of interest, a Multi-gap configuration does not seem to be mandatory for ion de- 
tection. A more practical RPC design could be probably based on 1-gap RPCs (single or mirrored 



[46]) 
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